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Abstract 

We show that localized iV-body soliton states exist for a quantum integrable derivative 
nonlinear Schrodinger model for several non-overlapping ranges (called bands) of the 
coupling constant rj. The number of such distinct bands is given by Euler's 0-function 
which appears in the context of number theory. The ranges of rj within each band can also 
be determined completely using concepts from number theory such as Farey sequences 
and continued fractions. We observe that iV-body soliton states appearing within each 
band can have both positive and negative momentum. Moreover, for all bands lying in 
the region r] > 0, soliton states with positive momentum have positive binding energy 
(called bound states), while the states with negative momentum have negative binding 
energy (anti-bound states). 
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1 Introduction 

Soliton states in integrable quantum field theory models in 1+1 dimensions have been 
studied extensively for many years [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. The quantum soliton 
states are usually constructed by using either the coordinate Bethe ansatz or the alge- 
braic Bethe ansatz. For an integrable nonrelativistic Hamiltonian, the coordinate Bethe 
ansatz can yield the exact eigenfunctions in the coordinate representation. If such an 
eigenfunction decays sufficiently fast when any of the particle coordinates tends towards 
infinity (keeping the center of mass coordinate fixed), we call such a localized square- 
integrable eigenfunction a quantum soliton state. It is also possible to construct quantum 
soliton states using the algebraic Bethe ansatz, by choosing appropriate distributions of 
the spectral parameter in the complex plane [2, 3, 4]. The stability of quantum soliton 
states, in the presence of small external perturbations, can be determined by calculating 
their binding energy. It is usually found that localized quantum soliton states of various 
integrable models, including the well known nonlinear Schrodinger model (NLS) and the 
sine-Gordon model, have positive binding energy [1, 2, 3, 4, 5]. 

In this paper, we will study the quantum soliton states of an integrable derivative 
nonlinear Schrodinger (DNLS) model [6, 7, 8, 9]. Classical and quantum versions of the 
DNLS model have found applications in different areas of physics like circularly polarized 
nonlinear Alfven waves in a plasma [11, 12], quantum properties of solitons in optical 
fibers [13], and some chiral Luttinger liquids which are obtained by dimensional reduction 
of a Chern-Simons model defined in two dimensions [14, 15]. It is known that the classical 
DNLS model can have solitons with momenta in only one direction, where the direction 
depends on the sign of the coupling constant; namely, the ratio of the momentum to the 
coupling constant is positive [16, 17, 18]. Here we want to investigate whether this chirality 
property of the classical solitons is preserved at the quantum level. The Lagrangian and 
Hamiltonian of the quantum DNLS model in its second quantized form are given by [6] 



where the subscripts t and x denote partial derivatives with respect to time and space 
respectively, i] (7^ 0) is the real coupling constant, and we have set the particle mass 
m = 1/2. The field operators ip(x, t), ip^x, t) obey the equal time commutation relations, 
[ijj(x,t),xjj(y,t)] = [i[)'(x, t), ip^(y, t)} = 0, and [ip(x, t), (y, 01 = h5{x — y). The quantum 
soliton states of this DNLS model have been constructed through both the algebraic Bethe 
ansatz [7, 8, 9] and the coordinate Bethe ansatz [6]. By applying the coordinate Bethe 
ansatz, it was found that quantum iV-body soliton (henceforth called iV-soliton) states 
exist for this DNLS model provided that r] lies in the range < \r]\ < tan (j^- Moreover 
it was observed that, similar to the classical case, such iV-soliton states can have only 
positive values of P/rj, where P is the momentum [6]. However, it was found recently 
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that soliton states can exist even with P/rj < provided that tan (j^j < \r)\ < tan [jv^ij 
[10]. These soliton states have the surprising feature that their binding energy is negative. 
This naturally leads to the question of whether there are additional ranges of values of i] 
for which there are iV-solitons in the quantum DNLS model, and if so, what the values of 
momentum and binding energy of those solitons are. 

In this paper, we re-investigate the ranges of values of rj for which localized quantum 
iV-soliton states exist in the DNLS model. In Sec. 2, we apply the coordinate Bethe ansatz 
and find the condition which the Bethe momenta have to satisfy in order that a quantum 
iV-soliton state should exist. In Sec. 3, we find that there are certain nonoverlapping 
ranges of rj, called bands, in which solitons can exist. The iV-solitons of DNLS model 
found earlier [6,10] are contained only within the lowest band. Thus the existence of 
solitons in higher bands is a novel feature of DNLS model which is revealed through 
our present investigation. We also apply the idea of Farey sequences in number theory 
to completely determine the ranges of all bands for which iV-soliton states exist for a 
given value of N. In Sec. 4, we show that the solitons appearing within each band can 
have both positive and negative values of P/rj] these have positive and negative binding 
energies respectively, and we call them bound and anti-bound states. In Sec. 5, we will 
use another concept from number theory, that of continued fractions, to give an explicit 
expression for the end points of the bands. We will also address the inverse problem of 
finding the values of N for which iV-soliton states exist for a given value of 77. We make 
some concluding remarks in Sec. 6. 



2 Conditions for quantum iV-soliton states in DNLS 
model 

To apply the coordinate Bethe ansatz, we separate the full bosonic Fock space associated 
with the Hamiltonian (1.1) into disjoint iV-particle subspaces \S N ). We want to solve the 
eigenvalue equation H\Sn) = E\Sn). The coordinate representation of this equation is 
given by 

H N T N (X!,X 2 , • • • ,x N ) = E r N (xi,x 2 , • ■ ■ ,x N ) , (2.1) 
where the iV-particle symmetric wave function tn(xi,x 2 , ■ ■ ■ ,xn) is defined as 

t n (x!,x 2 ,- ■ ■ ,x N ) = — L= (0|V>(zi) •• -ip(x N )\S N ) , (2.2) 



and H N , the projection of the second-quantized Hamiltonian H (1.1) on to the N particle 
sector, is given by 

H N = -h* £ — + 2ih* V E S( Xl -x m ) {q-+q-)- (2-3) 

j=l ux j l <m ux l ux m 
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It is evident that H N commutes with the total momentum operator in the iV-particle 
sector, which is defined as 

N Q 

p N = -in £ — . (2.4) 

j=l ux j 

Note that the Hamiltonian and momentum enjoy the scaling property Hn — > X 2 Hn 
and Pjy — > APjv if all the coordinates Xi — > Xj/A. Given any one eigenfunction of if at 
and Pjv, therefore, we can find a one-parameter family of eigenfunctions by scaling all 
the Xi. We also observe that if at remains invariant while P/v changes sign if we change 
the sign of 77 and transform all the Xi — > —Xi at the same time; let us call this the parity 
transformation. Hence it is sufficient to study the problem for one particular sign of rj, 
say, rj > 0. The eigenfunctions for rj < can then by obtained by changing Xj — > —xf, 
this leaves the energy invariant but reverses the momentum. 

Next, we divide the coordinate space R N = {x\,x 2 , ■ • • xn} into various iV-dimensional 
sectors defined through inequalities like x u (i) < x u (2) < •• • < x u (n), where u(l),u(2), 
■ ■ ■ ,u(N) represents a permutation of the integers 1, 2, • • • , N. Within each such sector, 
the interaction part of the Hamiltonian (2.3) is zero, and the resulting eigenfunction 
is just a superposition of the free particle wave functions. The coefficients associated 
with these free particle wave functions can be obtained from the interaction part of the 
Hamiltonian (2.3), which is nontrivial at the boundary of two adjacent sectors. It is 
known that all such necessary coefficients for the Bethe ansatz solution of a iV-particle 
system can be obtained by solving the corresponding two-particle problem [19]. Let us 
therefore construct the eigenfunctions of the Hamiltonian (2.3) for the two-particle case, 
without imposing any symmetry property on T2(x±, X2) under the exchange of the particle 
coordinates. For the region x\ < x 2 , we may take the eigenfunction to be 

t 2 (xi,x 2 ) = exp {i(kiXi + k 2 x 2 )} , (2.5) 

where k\ and k 2 are two distinct wave numbers. Using Eq. (2.1) for N = 2, we find that 
this two-particle wave function takes the following form in the region x\ > x 2 : 

t 2 (xi,x 2 ) = A(k 1 ,k 2 )exp {i(k±xi + k 2 x 2 )} + B(ki, k 2 ) exp {i(k 2 xi + k±x 2 )} , (2.6) 

where the 'matching coefficients' A(k± : k 2 ) and B(ki, k 2 ) are given by 

A(h, k 2 ) = h - h + lV [ kl + h) , B( kl , k 2 ) = 1 - A( kl , k 2 ) . (2.7) 

K\ — K 2 

By using these matching coefficients, we can construct completely symmetric iV-particle 
eigenfunctions for the Hamiltonian (2.3). In the region X\ < x 2 < ■•• < x^, these 
eigenfunctions are given by [6, 19] 

t n (x 1 ,x 2 ,---,x n ) = E II -irrrf^ P.fi)^),-,^)^!.^,---,^) , (2-8) 

u> \l< m A{K m) Ki) I 
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where 



P^(i)M^,-MN)(xi,x 2 ,- ■ ■ ,xn) = exp H h ^(ao^aO) • (2.9) 

In the expression (2.8), the k n 's are all distinct wave numbers, and XL implies summing 
over all permutations of the integers (1, 2, ....TV). The eigenvalues of the momentum (2.4) 
and Hamiltonian (2.3) operators, corresponding to the eigenfunctions tn(xi,X2, ■ ■ ■ , %n), 
are given by 

N 

P N t n (x 1 ,x 2 ,- ■ ■ ,x N ) = h(J2 k i) t n(xi,x 2 ,---,x n ) , (2.10a) 

3=1 
N 

H N t n (xi,x 2 ,---,x n ) = h 2 (^kfj t n (xi,x 2 ,- ■ ■ ,x N ) . (2.106) 

j'=i 

The wave function in (2.8) will represent a localized soliton state if it decays when 
any of the relative coordinates measuring the distance between a pair of particles tends 
towards infinity. To obtain the condition for the existence of such a localized soliton state, 
let us consider the following wave function in the region x\ < x 2 < • • • < x^ : 

N 

Pi,2,-,n( X!,x 2 , ■ ■ -,x N ) = exp (i^kjXj) . (2.11) 

3=1 

As before, the momentum eigenvalue corresponding to this wave function is given by 
ftX^Li kj. Since this must be a real quantity, we obtain the condition 

N 

Efc = °> ( 2 - 12 ) 

3=1 

where qj denotes the imaginary part of kj. The probability density corresponding to the 
wave function P\^,-,n{ Xi,x 2 , • • ■ ,x^ ) in (2.11) can be expressed as 

N-l r 

|Pi,2,...,jv( xi,x 2 ,---,x N )| 2 = exp { 2 (S*) ^ } ' ( 2 - 13 ) 

r=l j=l 

where the y r 's are the — 1 relative coordinates: y r = a; r+ i — x r , and we have used Eq. 
(2.12). It is evident that the probability density in (2.13) decays exponentially in the 
limit y r — > oo for one or more values of r, provided that all the following conditions are 
satisfied: 

N-l 

9i <0, gi + g 2 <0, , E ^<°- ( 2 - 14 ) 

3=1 

Note that the wave function (2.11) is obtained by taking u as the identity permutation in 
(2.9). However, the full wave function (2.8) also contains terms like (2.9) with oo represent- 
ing all possible nontrivial permutations. The conditions which ensure the decay of such 
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a term with a nontrivial permutation will, in general, contradict the conditions (2.14). 
Consequently, in order to have an overall decaying wave function (2.8), the coefficients of 
all terms Pu>(i),u>{2),-,u>(N)(xi, %2, • • • ,%n) with nontrivial permutations must be required to 
vanish. This leads to a set of relations like 

A(k u k 2 ) = 0, A(k 2: k 3 ) = 0, , A(k N _ u k N ) = 0. (2.15) 

Thus the simultaneous validity of the conditions (2.12), (2.14) and (2.15) ensures that 
the full wave function tn(xi,x 2 , • • • ,xn) (2.8) represents a localized soliton state. Using 
the conditions (2.12) and (2.15), we obtain an expression for the complex fc n 's in the form 

k n = X e-W 1 " 2 ")* , (2.16) 

where x is a re& l parameter, and is related to the coupling constant as 

= tan -1 (?7) . (2.17) 

To obtain an unique value of from the above equation, we restrict it to the fundamental 
region — | < 0(^ 0) < |. [Note that rj and have the same sign. Due to the parity 
symmetry mentioned above, we can restrict our attention to the range < < |]. 

In this context, it may be mentioned that a relation equivalent to (2.16) can also be 
obtained through the method of the algebraic Bethe ansatz, when quantum soliton states 
of DNLS model are considered [7, 8, 9]. 

Now, let us verify whether the fc n 's in (2.16) satisfy the conditions (2.14). Summing 
over the imaginary parts of these fc„'s, we can express the conditions (2.14) in the form 

^sm(70) sin r/ 7V _A^i > for l = x 2j ... N _ 1 / 218 n 
sin 

Thus, for some given values of 0, iV and x> a soliton state will exist when all the above 
inequalities are simultaneously satisfied. By using Eqs. (2.16) and (2.10), we obtain the 
momentum eigenvalue of such soliton state to be 

P = hx s Mm , (2 . 19) 

sm0 

and the energy to be 

h 2 X * sin(2iV0) 
E = sin(20) • (2 ' 20) 

The next section of our paper will be devoted to finding the ranges of values of where all 
the inequalities (2.18) are simultaneously satisfied for a given value of the particle number 
N. 
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3 Finding the values of where iV-soliton states exist 



In this section, we will study the values of where iV-soliton states exist for different 
values of N. For the simplest case N = 2, the condition (2.18) is satisfied when lies 
in the range < < | (— | < < 0) for the choice x > (x < 0). Thus any nonzero 
value of within its fundamental region can generate a 2-soliton state. Looking at the 
momentum P in (2.19), we see that the ratio P/0 > since and x have the same sign. 
Thus the chirality property of the classical solitons is preserved in the quantum theory 
for N = 2. 

We will now consider the more interesting case with N > 3. Due to the parity 
symmetry of the Hamiltonian in (2.3), we will henceforth assume that > 0. The 
inequalities in Eq. (2.18) can therefore be rewritten as 

X sin(/0) sm[(N-l)<f>] > for I — 1, 2, • • • ,N — 1 , (3.1) 

or 

X [ cos[(7V - 2l)<j>\ -- cos(iV0) ] > for I = 1, 2, • • • , N - 1 . (3.2) 
It is now convenient to consider the cases x > an d X < separately. 
For x > 0, Eq. (3.2) takes the form 

cos[(iV - 2Z)0] > cos(iV0) for I = 1, 2, • • • , N - 1 . (3.3) 

Let us now consider a value of of the form 

irn 

(PN,n = -Jjr , (3.4) 

where n is an integer satisfying 1 < n < N/2. Now, if n is odd, cos(iV07v,n) = —1 and it is 
possible that all the inequalities in (3.3) will be satisfied since —1 is the minimum possible 
value of the cosine function. However, a closer look reveals that all the inequalities in (3.3) 
are satisfied only if iV and n are relatively prime, i.e., if the greatest common divisor of N 
and n is 1. If N and n are not relatively prime, then let p be an integer (greater than 1) 
which divides both of them. Consider the integers N' = N/p and n' = n/p. Since n is odd, 
p must be odd and therefore n' is also odd. Similarly, N' = N/p is even (odd) if iV is even 
(odd), and therefore N — N' is an even integer which is less than N. We then see that there 
is an integer / for which N — 21 = N', and therefore cos[(A r — 2Z)0jv, n ] = cos(7rn') = — 1; 
hence that particular inequality in Eq. (3.3) will be violated. We therefore conclude that 
all the inequalities in (3.3) will be satisfied for = (p N>n and n odd, if and only if N and 
n are relatively prime. 

Similarly, we can consider the case x < 0- Eq. (3.2) then takes the form 

cos[(iV - 2/)0] < cos(iV0) for I — 1, 2, • • • ,N — 1 . (3.5) 

It is clear that all of these cannot be satisfied if iV is even; in particular the inequality with 
I = N/2 will fail. We therefore assume that iV is odd. Once again, we consider a value 
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of of the form given in (3.4) where n is now even. Now cos(iV07v,n) — 1 and there is a 
chance that all the inequalities in (3.5) will be satisfied since 1 is the maximum possible 
value of the cosine function. However, we again find that this is true only if N and n are 
relatively prime. If iV and n are not relatively prime, then let p be an integer (greater 
than 1) which divides both of them. Consider the integers N' = N/p and n! = n/p. Since 
iV is odd, p must be odd and therefore N' is also odd; hence N — N' is an even integer 
which is less than N. Similarly, n' = n/p is even since n is even. Thus there is an integer 
/ for which N — 21 = N', and then cos[(iV — 2Z)0jv, n ] = cos(7rn') = 1; hence that particular 
inequality is violated in Eq. (3.5). We therefore conclude that all the inequalities in (3.5) 
will be satisfied for = <pN,n and n even, if and only if N and n are relatively prime. 

Putting these statements together, we see that all the inequalities in (3.2) are satisfied 
for = <pN,n, if and only if N and n are relatively prime (with n odd for \ > 0, and 
n even for x < 0). By continuity, it then follows that all the inequalities will hold 
in a neighborhood of <pN,n extending from a value <f>N,n,- to a value 0jv,n,+, such that 
(f>N,n,- < 4>N,n < 4>N,n,+- We will call the region 

<t>N,n - < < <t>N,n,+ (3.6) 

as the band B^ >n . In this band, there is a soliton state with N particles. 

For a given value of N, the number of bands in which soliton states exists is equal 
to the number of integers n which are relatively prime to N and satisfy 1 < n < N/2. 
This is equal to half the number of integers which are relatively prime to iV and satisfy 
1 < n < N. The latter number is called Euler's 0-function or totient function <E>(iV) [20]. 
If iV has the prime factorization 

N = pTpT-'-pT , (3-7) 
where pi,P2) " • >Pn are a U prime numbers, then the totient function is given by 

k 1 

$(AQ = NH (1 - - ). (3.8) 

i=l Pi 

Thus the number of bands lying within the region rj > coincides with $(iV)/2. From 
the parity symmetry of Hamiltonian (2.3), it follows that the total number of bands lying 
within the full range of rj is given by $(iV). 

We now have to determine the end points <pN,n- and 4>N,n,+ of the band B^ n . The 
original inequalities in (3.1) show that the end points are given by of the form 

= y , (3.9) 

where j and / are integers satisfying 

1 < I < N , and j < l - (3.10) 
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(since <f) < n/2). In general, these j and I may not be relative prime numbers. However, 
one can always choose two relative primes j' and V satisfying the constraint (3.10) such 
that j'/l' = j/l. Thus the end points of the band B^,n are given by two rational numbers 
(f)/it of the form j/l (where j, / are relative primes satisfying the conditions in (3.10)) 
which lie closest to (and on either side of) the point 0jv,n/~ — n/N. The solution to this 
problem is well known in number theory and is described by the Farey sequences [20]. 

For a positive integer N, the Farey sequence is defined to be the set of all the fractions 
a/6 in increasing order such that (i) < a < b < N, and (ii) a and b are relatively prime. 
The Farey sequences for the first few integers are given by 

1 1 

12 1 

F 3 : °- i I I i 

1 3 2 3 1 

1112 3 1 
p ■ _______ 

1 4 3 2 3 4 1 

F 5 : ° i i i ? 1 3 2 3 4 1 
5 15435253451 K J 

These sequences enjoy several properties, of which we list the relevant ones below. 

(i) a/b < a'/b' are two successive fractions in a Farey sequence F N , if and only if 

a'b - ab' = 1 , and b , b' < N . (3.12) 
It then follows that both a and b' are relatively prime to a' and 6. 

(ii) For iV > 2, if n/N is a fraction appearing somewhere in the sequence F N (this implies 
that N and n are relatively prime according to the definition of i*jv), then the fractions 
ai/6i and 02/^2 appearing immediately to the left and to the right respectively of n/N 
satisfy 

ai , a 2 < n , and a x + a 2 = n , 

61 , b 2 < N , and b x + b 2 = N . (3.13) 

To return to our problem, we now see that the points <pN,n i n (3-4) (which lie in the 
bands Bn )H ) have a one-to-one correspondence with the fractions n/N, which appear on 
the left side of 1/2 within the sequence Fjv. Due to Eqs. (3.9) and (3.10), the end points 
of the band B N>n are given by 

<PN,n- — ~r~ , 
0\ 

<l>N,n, + = ^ , (3.14) 
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where a\jb\ and 02/62 are the fractions lying immediately to the left and right of n/N in 
the Farey sequence F^. These are the two unique fractions which lie closest to (and on 
either side) of n/N whose denominators are less than N (due to (3.13)). Therefore, as 
we move away from <p = (pN,n, one of the inequalities in (3.1) will be violated for the first 
time at these two points. (The property bi + b 2 = N shows that the same inequality in 
(3.1) is violated at the two end points of a given band B N;n ). 

The end points of a band given in Eq. (3.14) satisfy the property 

nb\ — Nai = 1 , 

nb 2 - Na 2 = - 1 , (3.15) 

due to the property (3.12) of Farey sequences. Given two integers N and n satisfying 
the conditions given above, one can numerically find the integers Oj and bi by searching 
for solutions of Eqs. (3.15) within the limits given in (3.13). For the case of lowest 
band (n = 1), these equations can also be solved analytically to yield: a± = 0, b\ = 
1, a 2 — 1, b 2 = N — 1. Thus, for every N > 3, the range of lowest band is given by 
< (p/n < 1/(N — 1) which agrees with the result obtained in Ref. [10]. It turns out that 
there is a way to analytically determine the integers and bi for the case of a general n 
by using the idea of continued fractions; this will be described in Sec. 5. 

In Table 1, we show the ranges of values of for which solitons exist for N = 2 to 9. 
In Fig. 1, we present the same information pictorially for iV going up to 20. 



N 


n 


Range of values 






of <p/n 


2 


1 


< 0/vr < 1/2 


3 


1 


< 0/tt < 1/2 


4 


1 


< 0/vr < 1/3 


5 


1 


< (f)/n < 1/4 


5 


2 


1/3 < 0/tt < 1/2 


6 


1 


< (f)/n < 1/5 


7 


1 


< 0/tt < 1/6 


7 


2 


1/4 < 0/vr < 1/3 


7 


3 


2/5 < <p/ir < 1/2 


8 


1 


< 0/vr < 1/7 


8 


3 


1/3 < 0/tt < 2/5 


9 


1 


< 0/vr < 1/8 


9 


2 


1/5 < 0/tt < 1/4 


9 


4 


3/7 < 0/vr < 1/2 



Table 1. The range of values of (p/n for which solitons exist for various values of N. 

Due to the relations in (3.15), we see that the width of the right side of the band B NjJl 
from (p N>n to (pN,n,+ is n/(Nb 2 ), while the width of the left side from (f>N,n - to 4>N >n is 
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7r/(iV&i). For later use, we note that each of these widths is strictly larger than tt/N 2 , 
since 61, b 2 < N. The total width Wn, u of the band B^^ is given by 

w "' = wi + m 2 ■ (3 ' 16 » 

For a given value of N, we will now find an expression for the total width of all the bands. 
Consider two different bands B Nn and B Nn i and their end points given by 



at 




n 




a 2 




< 




< 




h 


N 


















n' 




a' 2 




< 




< 




b[ 


N 


b' 2 



(3.17) 

where we know from (3.12) that n, n', b±, b 2 , b[, b' 2 are all relatively prime to N. We 
also know that n, n' < N/2; we will assume that n' < n. We will now show that b[ and 
b' 2 are not equal to either bi or b 2 . Eq. (3.12) implies that 

Ndi — nb\ = — 1 , 

Na 2 - nb 2 = 1 , 

Na[ - n'b[ = -1 . (3.18) 

We then see that 

N (ai — a[) = nb\ — n'b[ , 

N(a 2 + a' l ) = nb 2 + n'b[. (3.19) 

If b[ was equal to bi, we would get N(a\ — a[) = b'^n — n'), i.e., is a factor of b'^n — n'). 
Since b[ is relatively prime to N, this would mean that N must be a factor of n — n' [20]. 
But this is not possible since n — n' < N. Similarly, if b[ was equal to b 2 , we would get 
iV(a2 + a[) = b[(n + n'), i.e., is a factor of n + n' . But this is also not possible since 
n + n' < N. We therefore conclude that b[ is not equal to bi or b 2 . Similarly, we can show 
that b 2 is not equal to b± or b 2 . 

We thus have the result that for any value of N, the denominators of the end points 
of the various bands B^ n are all different from each other; we also know that they are 
all smaller than and relatively prime to N. Since the numbers of bands is equal to half 
the number of integers less than and relatively prime to N, and each band has two end 
points, the set of denominators of the end points of all the bands must be exactly the 
same as the set of integers less than and relatively prime to N. Eq. (3.16) now implies 
that the total width of all the bands is given by 

Wn = £ E j , (3.20) 
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where the sum runs over all values of I which are relatively prime to N and satisfy 
1 < I < N — 1. If iV is prime, we obtain 

W n = £ £ 7 - ^ [ ^-1) + 7] as N ^ oo , (3.21) 

where 7 = 0.57721 • • • is Euler's constant. We have numerically studied the behaviors of 
Wn and the sum I N = X^ =2 ^v for values of N up to 10 4 . We find that although Wn is 
a non-monotonic function of N, on the average it grows with N as In N/N; consequently, 
In grows as (In TV) 2 as N becomes large. [In the range N = 10 2 to 10 4 , we find that 
NW N /[n(ln(N - 1) + 7)] fluctuates between 0.255 and 1.001, while W N /(\nN) 2 = 1.179 
for N = 10 4 ]. 

We thus see that the allowed range of values of for which solitons exist goes to zero 
generically as \nN/N as iV — > 00. In contrast to this, the lowest band runs from to 
n/(N — 1) and therefore has a width of ix/(N — 1). This implies that the width of the 
lowest band becomes an insignificant fraction of the total width as N becomes large. 



4 Momentum and binding energy of a iV-soliton state 

In this section, we will calculate the momentum and binding energy for the iV-soliton 
states described above. We first look at the momentum of the solitons in a particular band 
B N , n using Eq.(2.19). The form of the end points given in (3.14) shows that sin(iV0) = 
at only one point in the band Bn,u, namely, at = (pN,n- In the right part of the band 
(i.e., from (f) N , n to 4>N,n,+)i the sign of sin(A r 0) is (— l) n . In the left part of the band 
(i.e., from <pN,n - to <f>N,n), the sign of sin(iV0) is (— l) n+1 . Now, the analysis given above 
showed that x has the same sign as (— l) n+1 . Hence the momentum given in (2.19) is 
positive in the left part of the band, negative in the right part of the band, and zero at 

= 4>N,n- 

Next, we look at the energy using Eq. (2.20). To calculate the binding energy, we 
consider a reference state in which the momentum P of the iV-soliton state given in (2.19) 
is equally distributed among the iV single-particle scattering states. The real wave number 
associated with each of these single-particle states is denoted by ko- From Eqs. (2.10a) 
and (2.19), we obtain 

= xpm, (4 .D 

TV sin 

Using Eq. (2.10b), we can calculate the total energy for the iV single-particle scattering 
states as 

E - h 2 Nk 2 _ ^Vsin 2 (iV0) 

^ ~ ni ™° ~ N sin 2 • [ ^ Z) 
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Subtracting E in (2.20) from E s in (4.2), we obtain the binding energy of the iV-soliton 

state as 

E M ,N) S E.-E = "V^WIjEWl _ . (4.3) 

sin0 L iVsin0 cos0 ' 

It may be noted that the above expression of binding energy remains invariant under the 
transformation — > — 0. 

Substituting iV = 2 in Eq. (4.3), we obtain E B (<j>,2) = 2^\ 2 sin 2 0. Thus we get 
E B (4>,2) > for any nonzero value of 0. Let us now look at the case N > 3. We can 
rewrite Eq. (4.3) in the form 

Eb[4 ,, N) = ^sinW) 
JM sin cos 

/(0,7V) = x [ sin(iV0)cos0 - iVcos(iV0) sin0 ] . (4.4) 

We will now prove that the function /(0, N) is positive in all the bands B N>n for all values 
of N and n. To show this, we add up all the inequalities given in (3.2), and use the 
identities 

N t cos K iv - m = , 

JT± sin 



sin[(iV - 1)0] = sin(iV0) cos - cos(iV0) sin . (4.5) 

We then find that /(0, N) > 0. Hence, E B given in (4.4) has the same sign as x sin(iV0). 
Following arguments similar to that of the momentum, we find that the binding energy 
is positive in the left part of each band (called bound states), negative in the right part 
(called anti-bound states), and zero at the point = <pN,n- 

We thus see that for > 0, the momentum and the binding energy are both positive 
in the left part of each band, and they are both negative in the right part. This is to 
be contrasted with the classical solitons in the classical DNLS model which always have 
positive momentum. If < 0, we can similarly show that solitons with positive values of 
Pj have positive binding energy, and solitons with negative values of P/ have negative 
binding energy. 

In Fig. 2, we show the binding energy E B in (4.3) as a function of 0/7T for three 
different values of N. (We have set h 2 x 2 = 1 in that figure). We see that E B is indeed 
positive (negative) in the left (right) part of each band. 



5 Continued fractions 



In this section, we will apply the technique of continued fractions to our problem. This 
will lead us to an explicit expression for the end points of the band B Nn for any value 
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of TV and n. We will also address the problem of determining the values of N for which 
iV-soliton states exist for a given value of <fi. It will turn out that continued fractions 
provide a way of finding some (but not necessarily all) values of N for which solitons exist 
for a given 0. 

We first discuss the idea of a continued fraction [20]. Any positive real number x has 
a simple continued fraction expansion of the form 

1 . . 

x = n H i — , (5.1) 

^1 ~~T , l 

where the n^s are integers satisfying n > 0, and rij > 1 for i > 1. The expansion ends 
at a finite stage with a last integer n k if x is rational. In that case, we can assume that 
the last integer satisfies n k > 2. (If n& is equal to 1, we can stop at the previous stage 
and increase rik-i by 1). With this convention for n k , the continued fraction expansion 
is unique for any rational number x. If x is an irrational number, the continued fraction 
expansion does not end, and it is unique. 

We will be concerned below with the continued fraction expansion for x = <fi/ix which 
lies between and 1/2; hence, no = 0. We will use the notation x = [0, ni, ri2, ri3, • • •] for 
such an expansion. 

Given a number x, the integers rii can be found as follows. We define x = x. Then 
n o — [ x o], where [y] denotes the integer part of a non-negative number y. We then 
recursively define = — rij), and obtain n i+ i = [x i+ i] for % = 0, 1,2, • • •. If we 

stop at the k th stage, we obtain a rational number r k = [n ,ni,n 2 , ■ ■ • ,n k ] which is an 
approximation to the number x. If we write r k = Pk/qk, where pk and q k are relatively 
prime, then we have the following three properties for all values of k > 1 [20] , 

<?fc < Qk+l , 
Pk+iQk ~ PkQk+i = (-1)^ , 

I a: - — I < \ . (5.2) 

Using these properties, we can now find the end points <f>N,n,± of the band B^, n . In Sec. 
3, we saw that this is equivalent to finding the fractions a\jb\ and 02/62 which lie to the 
immediate left and right of the fraction n/N in the Farey sequence F^. Let us suppose 
that n/N has a continued fraction expansion given by 

n 

— = [0,ni,n 2 , • • • ,n k -i,n k \ 

= [0,ni,n 2 ,---,n fc _i,n fe -l,l] , (5.3) 

where the expression in the second equation can be seen to be equivalent to the one in 
the first equation. Now consider the two fractions 

a 

- = [0,7ii,n 2 , • • • ,n fc _ij , 

Q 

- = [0,ni,n 2 , • • • ,n k -i,n k - 1] . (5.4) 
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On comparing these expressions to those given in (5.3) and using the properties of con- 
tinued fractions in (5.2), we see that 



b, d < N , 
nb — Na = (-if -1 , 

nd — Nc = {-l) k . (5.5) 

We now see from Eq. (3.12) that if k is odd, then a/b and c/d are to the immediate left 
and right respectively of n/N in the Farey sequence F N , while if k is even, then a/b and 
c/d are to the immediate right and left respectively of n/N in the sequence F N . Hence 
the end points (pN^-/ 71 and <f> Njni+ /n are given by a/b and c/d respectively (in Eq. (5.4)) 
if k is odd, and vice versa if k is even. Thus the continued fraction expansion gives a 
convenient way of determining <pN,n,± if N and n are given. 

As an explicit example, consider the value 0/7T = 3/8. The continued fractions ex- 
pansion of this is given by 4>/it = [0,2,1,2] = [0,2,1,1,1]. Eq. (5.4) then gives the 
two neighboring fractions in the Farey sequence F 8 to be a/b = [0,2,1] = 1/3 and 
c/d = [0,2,1,1] = 2/5. The band B 8 , 3 therefore lies in the range 1/3 < 0/tt < 2/5 
as shown in Table 1. 

Next, we study the question of determining the values of N for which TV-body solitons 
exist if a value of is given. Suppose that we know the continued fraction expansion 

- = [0,ni,n 2 , • • •] . (5.6) 

7T 

Suppose that we stop at any any point in this expansion, say, at the k th stage, and we get 

— = [0,rai,ra2, • • • ,n k ] ■ (5.7) 



Then we know from (5.2) that 



- - — I < 4 ■ (5-8) 



We now recall from Sec. 3 that both the right and the left part of the band B qk Pk have 
widths which are larger than \jq\. Eq. (5.8) therefore implies that (fi/n must lie within 
the band B qktPk . We have thus found a value of N = q k for which we know that a iV-body 
soliton state exists for the given value of 0. We can generate several such values of iV by 
stopping at different stages k in the expansion given in (5.6). 

If 0/7r is rational, the continued fraction expansion stops at a finite stage, so we only 
obtain a finite number of values of N in this way. This can also be seen directly from Eq. 
(3.1). If <f)/n = p/q is rational, then at least one of the inequalities in (3.1) will be violated 
if iV > q. Moreover, one gets k n = k n+q from Eq.(2.16) by putting (p/n = p/q. Thus at 
least two wave numbers would coincide when N > q. As a result, the eigenfunction (2.8) 
becomes trivial in this case. We thus conclude that if (p/n is rational, there is only a finite 
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number of values of N for which a iV-body soliton state exists. On the other hand, if 
0/7T is irrational, then the expansion in (5.6) does not end, and we can use the procedure 
described above to find an infinite number of possible values of N for which a iV-body 
soliton exists. 

Given a value of 0, the procedure described above will find values of N for which 
we can prove that a iV-body soliton exists. This does not rule out the possibility that 
there may be other values of N for which such solitons exits. [To give a simple example, 
consider the value = 7r/10. Using Eq. (3.1), we see that there are bound states for all 
values of N from 2 to 10. However, the continued fraction technique only gives the value 
of N = 10 for which a soliton exists]. The continued fraction technique therefore provides 
a sufficient but not necessary condition for finding the desired values of N. 



6 Conclusion 

By applying the coordinate Bethe ansatz, we have investigated the range of the coupling 
constant r\ for which localized quantum iV-soliton states exist in the DNLS model. Using 
the ideas of Farey sequences and continued fractions, we have given explicit expressions 
for all the allowed ranges (bands) of i] for which iV-body solitons exist. The continued 
fraction also provides a way of finding some values of N for which soliton states exist for 
a given value of rj. If = tan _1 (?7) is rational (irrational), the number of values of N for 
which soliton states exist is finite (infinite). 

We find that for N > 3, the iV-body solitons can have both positive and negative 
momentum. Thus the chirality property of classical DNLS solitons is generally not pre- 
served at the quantum level. We also calculated the binding energy for the soliton states. 
We find that solitons with positive values of P/rj form bound states with positive binding 
energy. Solitons with negative values of P/rj have negative binding energy and hence 
form anti-bound states. The anti-bound solitons would be expected to be unstable in the 
presence of external perturbations. As a future study, it may be interesting to calculate 
the decay rate of such quantum solitons by introducing small perturbations in the DNLS 
Hamiltonian. 

In the continuum version of the DNLS model (this can be obtained from the iV-body 
Schrodinger model by taking the limit N — > oo), it is known that solitons exist only if 
< N\(j)\ < 7C [6, 18]; this corresponds to the lowest band. One might ask whether there 
is a continuum version of the solitons in the higher bands which we have found in this 
paper. This may be an interesting problem for future study. It is possible that only the 
lowest band has a counterpart in the continuum theory. Note that if lies in the lowest 
band, then the phases of k n in Eq. (2.16) vary slowly with n; the spacing between two 
neighboring phases is of order 1/N. This may be necessary to ensure that the continuum 
limit exists. For the generic band, where is of order 1 rather than of order 1/N, the 
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phases of neighboring values of k n vary rapidly; this may make it difficult to take the 
continuum limit. 

To conclude, we have shown that the soliton structure in the DNLS model is much 
richer than in the well known NLS model. In the latter model, the interaction part of the 
Hamiltonian is given by nJ2i< m H x i ~ x m) in the coordinate representation. This model 
has iV-soliton states for any value of N > 2 and any value of pL < 0; all the solitons have 
positive binding energy. In contrast to this, the DNLS model allows soliton states in only 
certain bands of the coupling constant for N > 3, and these bands form intricate number 
theoretic patterns. 
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Figure 1: The values of 4>/tt for which iV-body soliton states exist for various values of N. 
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Figure 2: The binding energy E B of the iV-body soliton states as a function of 4>/ti 
for three different values of N . Positive and negative values of E B denote bound and 
anti-bound states respectively 
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